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1 INTRODUCTION 

• I— I , 

. Initiated with the work on atmospheric convection by Ogura 
5_j ■ & Phillips (1962) and Gough (1969), the anelastic approx- 

I imation has been widely used in astrophysics to describe 
compressible convection. For example, it has been first ap- 
plied by Latour et al. (1976) and Toomre et al. (1976) to 
investigate the non-rotating and non-magnetized convective 
envelopes. Then, Oilman & Glatzmaier (1981) and Glatz- 
maier & Gilman (1981a, b,c) extended these studies to the 
rotating case whereas the MHD case has been considered 
by Glatzmaier (1984,1985a,b), Lantz & Sudan (1995) and 
Lantz & Fan (1999). 

In all these studies, this approximation has been pre- 
ferred to an other famous one in fluid dynamics, the so-called 
Boussinesq approximation where the velocity is assumed to 
be solenoidal, i.e. V - v = (Spiegel & Veronis 1960). As the 
density varies on several orders of magnitude in a convection 
zone, using the anelastic approximation is indeed more rele- 
vant and leads to V-(pu) = 0. However, both approximations 
filter out short-period acoustic oscillations by assuming that 
the Mach number of the convection is small so that 'the task 
of numerical solution will not be complicated by the need to 
resolve very rapid time variations' (Latour et al. 1976). 

More recently, the anelastic approximation has been 
used in the field of stellar oscillations. Waves can propagate 
over long distances in star interiors and large variations of 
density need also to be taken into account. Thus, as in stel- 
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lar convection studies, the Boussinesq approximation is too 
restrictive and not adapted for this problem. However, since 
acoustic waves are filtered out, only the low- frequency part 
of the oscillations spectrum is captured. As this part of the 
spectrum also corresponds to the most perturbed one when 
rotation acts, this approximation is very attractive when 
low-frequency modes of a rotating star are considered. 

Understanding the spectrum of rotating stars is indeed 
a difficult problem when the rotation period is of the same 
order than the oscillation one. In this case, the usual per- 
turbative theory of Ledoux (1951) reaches its limits and the 
eigenvalue problem needs to be solved non-perturbatively 
(Dintrans & Rieutord 2000, hereafter referred to as Paper 
I). But rotation introduces new computational challenges: 
(i) the strong Coriolis coupling between normal modes of 
different degrees {£,£ ± 1) leads to a large system of cou- 
pled differential equations; (ii) some singularities due to the 
presence of wave attractors emerge in the region a ~ 20,, a 
being the wave frequency and fl the uniform rotation rate 
of the star. Applying the anelastic approximation does not 
remove these singularities but decreases the size of the nu- 
merical problem since acoustic quantities disappear (see Pa- 
per I). We note that Berthomieu et al. (1978) used the same 
arguments when considering a hierarchy between the dis- 
placement and pressure components, this hierarchy leading 
to the same result as the direct use of the anelastic approx- 
imation. 

Contrary to stars, the mass of the Earth is not strongly 
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concentrated near its centre. Only 30% of the whole mass 
is indeed inside the half radius whereas this ratio is around 
90% for the Sun (see e.g. Stix 1989) . It means that the vari- 
ations in density in the Earth's outer fluid core generated 
by the propagation of seismic P-waves may produce non 
negligible variations in the gravitational potential (f) and ex- 
plains why geophysicists do not wish to use the Cowling 
approximation in their description of P-waves. To take into 
account this effect but to simplify the equations for low fre- 
quency modes Smylie & Rochester (1981), Smylie, Szeto & 
Rochester (1984) introduced the 'subseismic' approximation 
where self-gravity is included. 

Recently, De Boek, Van Hoolst & Smeyers (1992) ap- 
plied this approximation to describe the low-frequency g- 
modes of non-rotating stars. Starting from the same equa- 
tions as Smylie & Rochester (1981), they found an analytic 
expression for the asymptotic low frequencies which slightly 
differs from the standard one deduced by Tassoul (1980) 
in the Cowling approximation. By applying this expression 
to the low-frequency oscillations of both homogeneous and 
polytropic (n — 3) star models, De Boek et al. (1992) con- 
cluded that the subseismic approximation is more accurate 
than the standard one in the polytropic case but less accu- 
rate in the homogeneous case. 

However, as will be shown below, the subseismic ap- 
proximation is not a consistent approximation of the equa- 
tions of motion; we will show indeed that the proper way 
to simplify the equations in order to deal with the low fre- 
quency modes is to use the anelastic approximation and that 
this approximation also implies Cowling's approximation. 

Hence, after establishing the complete set of equations 
describing linear oscillations of a stratified compressible fiuid 
and the appropriate boundary conditions (Section 2), we 
focus on the similarities and differences of the anelastic 
and subseismic approximations (Section 3); in particular, we 
show that both imply Cowling's approximation and demon- 
strate the inconsistency of the subseismic approximation. 
Then we compare the two approximations on the test cases 
used by De Boek et al. (1992), namely, the homogeneous 
and polytropic star models (Section 4). Comparison of the 
resulting approximate eigenfrequencies with the exact ones 
(computed numerically for the polytrope but analytically 
for the homogeneous model) shows clearly the better effi- 
ciency of the anelastic approximation. Finally, we conclude 
in Section 5 with some outlooks of our results. 



2 THE BASIC EQUATIONS 

Assuming a time-dependence of the form exp{iat), the gov- 
erning equations that describe the adiabatic oscillations of 
a spherically non-rotating star are given by 



where P',p' and (p' respectively denote the Eulerian fluc- 
tuations of pressure, density and gravitational potential 
whereas SP, 5p and ^ are the Lagrangian fluctuations of pres- 
sure, density and displacement such as 



5P = P' 



dP . 
dr 



5p- 



I .dp 



with a pressure gradient satisfying the hydrostatic equilib- 
rium dP/dr — —pg. Also, p is the equilibrium density, g the 
gravity and 7 = (9 In P/d In p) s the flrst adiabatic exponent . 
Finally, c? and A'^^ respectively denote the squares of sound 
speed velocity and Brunt- Vaisala frequency such as 



2 P 
c =7 — , 

P 



^2 / 1 dlnP dlnp 
^ 7 dr dr 



(5) 



In order to obtain a well-posed eigenvalue problem with an 
eigenvalue cr^, boundary conditions should be stipulated at 
the star centre and surface. At the centre, we impose the 
regularity of all pulsation quantities such as (see e.g. Unno 
et al. 1989) 



4r oc r , P cc r , 



oc r , 



(6) 



where £ denotes the spherical harmonics degree of the eigen- 
mode. At the surface, which is assumed to be free, we im- 
pose the following condition coupling and </!>' (Ledoux and 
Wahaven 1958) 



dr I + T*^'^^^ ^ -4^G(pC.)«. 

R 



(7) 



Finally, we adopt the classical mechanical outer condition 
for the pressure which reads 



5P = at r ^ R. 



(8) 



3 THE SUBSEISMIC AND ANELASTIC 

APPROXIMATIONS APPLIED TO STELLAR 
OSCILLATIONS 

3.1 The common properties 

Since the work of Cowling (1941), it is well known that g- 
modes are mainly transverse whereas p-modes are predom- 
inantly radial. The Eulerian pressure and density perturba- 
tions are thus small for a g-mode and both subseismic and 
anelastic approximations take advantage of that by assum- 
ing that the Lagrangian pressure fluctuation SP is only due 
to the radial displacement one; i.e. the Eulerian fluctuation 
P' does not contribute and we have 

"^^^r = -Pg^r- (9) 



5P = p' + 

dr 



dr 



In this case, Eqs. (^), (^ and (^) read now 



p' + V • (pC) = 0, 



a'f = V ( — + - —SPer 



SP = cSp, 
V^0' = 47rG'p', 



pg 



(1) 

(2) 
(3) 

(4) 



P = P^r, 

9 



N 



(10) 

(11) 
(12) 



At this stage, two important consequences appear: 



The anelastic and subseismic approximations for low-frequency g-modes 3 



(1) As observed by De Boek et al. (1992), the radial La- 
grangian displacement necessarily vanishes at the surface. 
To show this, we take the radial component of Eq. ( p^ 



The problem arises from the fact that the Brunt- Vaisala fre- 
quency diverges as 1/c^ at the star surface (see Eq. |5|). Thus, 
as diverges to infinity at r = i?, must vanish at this 
point to avoid a singular radial displacement. In the com- 
plete case (i.e. when acoustic waves are included), we note 
that it is the mechanical condition 5P — which permits 
a finite radial displacement at the surface, as shown by the 
radial component of Eq. (b|) 



2f 



d_f P' 

dr 



+ > 



pg 



SP, 



where the term N"^ /{pg) diverges as l/{p(?). Thus, the sub- 
seismic and anelastic ^, -eigenfunctions always have one node 
less than their corresponding unapproximated eigenfunc- 
tions for which £,r{R) is finite. 

(2) The second consequence, which is not straightfor- 
ward, is that the use of the subseismic or anelastic approx- 
imation necessarily implies Cowling's approximation where 
the perturbations (j)' are neglected (Cowling 1941). To show 
this, we first take the curl of Eq. (lid) and obtain 

a^V X e = V X (iV^Crfir). 

Furthermore, combining Eqs. (^ and ( pl| ) allows us to find 
the following subseismic form for the equation of mass con- 
servation 



— P^r + V ■ (pO = 



(13) 



whereas, as shown in the introduction, the anelastic coun- 
terpart of this equation reads 



V • (pO = 



dr 



(14) 



The two different forms (|l^) and (^^ may be formally writ- 
ten £{£,) — where £ denotes a differential operator de- 
pending on the approximation we use. We have thus to deal 
with the new following system 



O-V X 5 = V X (iV^Cr 



= 0, 



V^0' = 4ttG 



(15) 



(16) 



(17) 



whereas the surface boundary conditions now read 



/ 4- 1 

+ -Tr-'t>'iR) = 0. 



C,.(i?)=0 and 



R 



It is thus clear that the Poisson equation ( p^ decouples 
from Eqs . ( p^ ) and (|l^). In other words, the knowledge of 
Eqs. disl ) and (^6|) with the boundary condition ^r{R) = 
is sufficient to find the set of eigenfrequencies a^. The grav- 
itational potential perturbation </>' does not act on the de- 
termination of and can be seen as 'slaved' to the via 



Eq. (0)0 We note that Smylie & Rochester (1981) do not 
clearly emphasize this point and just noticed that the <j>'- 
perturbations can be determined after solving their subseis- 
mic equations for velocity and reduced pressure. 

3.2 The main difference 

As shown above, the subseismic and anelastic approxima- 
tions do not lead to the same equation of mass conservation 
since the term stemming from the density perturbations is 
kept in the subseismic case (see Eqs. |l| and |l^. We will 
now show that the subseismic form of this equation is in 
fact not correct because it is not consistent with the basic 
assumption from Eq. (^. 

First of all, we recall that both approximations assume 
that the Eulerian pressure perturbation does not contribute 
to the Lagrangian one, that is 



P' ^ dlnP 



(18) 



where Hp = {—din P/dr)~^ denotes the pressure scale 
height. Moreover, we can relate the Eulerian pressure and 
density perturbations using Eq. (Bl) as 



p 



9 



1 



1 



where Hp = (— dlnp/dr)^^ denotes the density scale height 
and where Eq. has also been used. Therefore the basic 
assumption from Eq. ( p^ leads to 



P_ 
P 



^ _ ^ u « ^ 

H, jHpj'^^^jHp 



p tip 



This last equation is important since it shows that we can 
also neglect the Eulerian contribution p' in the Lagrangian 
perturbation Sp when Eq. (^) is assumed, i.e. we have 

5p^p' + ^/jL^^/^. 

dr dr 

In other words, if we just take into account the contribution 
stemming from the equilibrium pressure gradient in SP, we 
should also do the same with the contribution coming from 
the equilibrium density gradient in 5p. Hence 

(5P ~ fr— ;— is equivalent to 5p~ (20) 
dr dr 

We can now easily understand why the subseismic equa- 
tion of mass conservation is not correct. Eq. rtlh indeed gives 



p' + V ■ {pi) = 



<5p-hpV-C = 0. 



In order to be consistent with the basic assumption from 
Eq. (|20[), the equation of mass conservation now reads 



dr 



dlnp 
dr 



and we recover the good anelastic form (|l^. 

The subseismic form of mass conservation ( [l^ of Smylie 
& Rochester (1981) is therefore clearly inconsistent because 

* It explains why Smeyers et al. (1995), starting from the com- 
plete case, found the same asymptotic development of as those 
derived by Tassoul (1980) under the Cowling approximation. 
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of its incompatibility with the hypothesis (|9|). As a conse- 
quence, the subseismic approximation is necessarily less ac- 
curate than the anelastic one as we will now show. 



4 RESULTS 

In order to test the efficiency of both approximations, we 
applied the following procedure to the oscillations of the 
homogeneous compressible star (i.e. a polytrope of order 0) 
and the polytrope n = 3: 

(1) we first calculated the exact eigenfrequencies of the 
full case, that is we solved the set of equations (1-4) using 
boundary conditions (6-8). This has been achieved analyti- 
cally for the polytrope n = (Pekeris 1938) and numerically 
for the polytrope n = 3. 

(2) we next calculated the subseismic and anelastic eigen- 
frequencies by solving the system consisting in Eq. (|l^) , the 
equation of mass conservation ( p^ or (^^ and finally the 
boundary condition ^r{R) ~ discussed in Section 3.1. Once 
again, analytical expressions have been used in the homoge- 
neous case whereas the polytrope n = 3 has been computed 
numerically. The efficiency of both approximations has been 
deduced by a direct comparison between these approximated 
eigenfrequencies and the previous exact ones. 



4.1 The homogeneous polytrope 

As a first example, we concentrate on the unstable second- 
degree g~ -modes of the polytrope n — 0. Pekeris (1938) first 
studied the nonradial oscillations of this model and found 
the exact g~ -eigenfrequencies as 



Table 1. Squares of exact ui'^ 



uj' = Att - v^Af^ + A, 
where A = ^(^ -f 1) and 



(fc = 0,l,2,...). 



- 2. 



Here lu denote the dimensionless eigenfrequencies, which for 
a star of mass M and radius R, are related to the cr^-ones 
by 

2 GM 2 

Using an equivalent method as the one used by Pekeris 
(i.e. the search of a power-serie solution), we calculated the 
analytic expressions of the subseismic and anelastic eigen- 
frequencies of the homogeneous polytrope and found (see 
appendix A) 



2 



2 _ 2 



7 A - (£ + 2k') {£ + 2k' + 1 + ^)'' 
A 



7A-(£ + 2fc')(^ + 2fc'-f-l)' 



(21) 



(22) 



where (fc' = 1,2, . . .). 

We summarized in Table 1 the exact eigenfrequencies 
calculated for the homogeneous star with 7 = 5/3, 1 = 2 and 
various orders fc, whereas Fig. |l| gives the relative errors (in 
percents) made with both approximations for the first thirty 
g~ -modes. The anelastic approximation clearly appears to 



subseismic ui'^ , and anelastic 

sub 



eigenfrequencies of the homogeneous polytrope. 







a;2 , X 10^ 

sub 


u)'^ , X 10^ 

and 




-7.251703 
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-6.923077 


Ss 


-3.613671 
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-1.212121 


sFs 
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-0.608108 


g20 


-0.365629 
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-0.243308 
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Figure 1. Rates of relative errors of the subseismic (Eq. |2l| ) and 
anelastic (Eq. [2^ eigenfrequencies as the order k increases for the 
homogeneous star model. 

be more precise than the subseismic one since, for example, 
the anelastic eigenfrequency at fc = 30 is about twenty times 
more precise than the subseismic one. 

4.2 The polytrope n = 3 

4-2.1 Numerics 

We next calculated the second-degree g+-mode of the poly- 
trope n = 3. As already mentioned, analytic expressions of 
eigenfrequencies are not known for this model then we com- 
puted numerically the eigenfrequencies of g"''-modes with 
the complete set of equations and their approximated sub- 
seismic and anelastic counterparts. As in our preceding work 
(Dintrans, Rieutord & Valdettaro 1999; paper I), we used an 
iterative eigenvalue solver based on the incomplete Arnoldi- 
Chebyshev algorithm. This numerical method differs from 
the classical ones where eigenfrequencies are determined by a 
direct integration of the governing equations using either re- 
laxation methods (Osaki 1975) or shooting methods (Hansen 
& Kawaler 1994) . Our eigenvalue formulation is presented in 
appendix B, with a discussion on the tricky problem of the 
degeneracy of the eigenvalue equations at the star surface 
when the subseismic or anelastic approximation are used. 

Eigenfrequencies are solutions of a generalized eigen- 
value problem 
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Table 2. Squares of the true A^, subseismic A^^^^ 
cigcnfrequencies of the polytrope n = 3. 



and anelastic 



and 





A2 X 10^ 


A2 , X 10^ 

sub 


A2 , X 10^ 

and 


g3 


1.121064 


1.269076 


1.210566 


g5 


0.576059 


0.627794 


0.603954 


gio 


0.198396 


0.208345 


0.202646 


gl5 


0.099755 


0.103233 


0.100973 


g20 


0.059984 


0.061587 


0.060453 


g25 


0.040039 


0.040911 


0.040259 


g30 


0.028631 


0.029154 


0.028744 



where Ma,M b denote two difFerential operators discretized 
on a Gauss- Lobatto grid associated with Chebyshev polyno- 
mials and H is the eigenvector associated with A^. We also 
note that the regular singularities appearing in the subseis- 
mic or anelastic equations at the surface are replaced by the 
boundary conditions (see appendix B.2). 



7 0. 



■< 1.0 




0. 1 



5 10 15 20 25 30 



Figure 2. Rates of relative errors of computed eigenfrequencies 
for the polytrope n = 3. 



4.2.2 Results 

We summarized in Table 2 the computed eigenfrequencies 
(in units of A-nGpc) for the polytrope n = 3 with 7 = 5/3, 
1 = 2 and various radial order^. As in the homogeneous 
case, the anelastic approximation appears to be superior to 
the subseismic one. This is illustrated by Fig. ^ where we 
show that the anelastic relative errors are always smaller 
than the subseismic ones. Still for k = 30, the agreement 
with the true eigenvalue is of about five times better with 
the anelastic approximation than with the subseismic one. 

We note, however, that this accuracy difference tends 
to be less pronounced than for the homogeneous model. In 
fact, the comparison of Fig. |l] and ^ shows that the sub- 
seismic relative errors are the same in both cases whereas 
the anelastic ones increase in the polytropic case. Therefore 
dealing with a model with important density variations at- 
tenuates the differences between the two approximations. 



5 CONCLUSION 

We studied the properties of the subseismic and anelastic 
approximations devised for low-frequency g-modes of stars. 
These approximations both assume that the Eulerian per- 
turbations of pressure do not contribute to their Lagrangian 
parts. As a consequence, acoustic waves are filtered out. 
However, the equation of mass conservation differs since 
the subseismic approximation keeps the contribution com- 
ing from the Eulerian density fluctuation. We showed that 
this is incorrect because it is inconsistent with the neglect 
of the Eulerian pressure perturbation. 

As an illustration we first applied both approximations 
to the low- frequency g~ -modes of the homogeneous star. In 
this case, we found an analytic expression for the subseismic 

t To obtain eigenfrequencies in the same units than in the homo- 
geneous case, we should multiply A^ by a;i/g ~ 162.547; i.e. we 
have up' = (xi/g)A^. 



and anelastic eigenfrequencies using power-series solutions. 
We thus showed that the anelastic approximation is more 
precise than the subseismic one as expected. 

As a second test, we considered the low- frequency g"''- 
modes of the polytrope n = 3. No analytic solutions can 
be found for this model and eigenfrequencies have been 
computed numerically using an iterative solver. We thus 
computed the subseismic and anelastic eigenfrequencies and 
compared them to the exact ones. As in the homogeneous 
case, the anelastic approximation is more accurate than the 
subseismic one. 

We also showed that Cowling's approximation is a nec- 
essary consequence of the anelastic approximation; it is 
therefore unnecessary to enforce this approximation since 
perturbations of self-gravity will automically decouple from 
the set of equations and will not influence the eigenfrequen- 
cies. 

As already mentioned in the introduction, an obvious 
application of the anelastic approximation concerns the low- 
frequency oscillations of rapidly rotating stars. Applying this 
approximation indeed removes the acoustic quantities and 
dramatically decreases the size of the numerical problem. An 
example is given in Paper I with the study of the oscillations 
of a rapidly rotating 7 Doradus-type star, for which periods 
of oscillations are of the same order as that of rotation (i.e. 
around one day). 
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APPENDIX A: ANALYTIC EXPRESSIONS FOR 
THE HOMOGENEOUS STAR MODEL 



Al The subseismic case 
We start from Eqs. (|l|) and (jl|) 
cr^V X f = V X (iV^Cre,.), 

where g,c? and A''^ are given by 



(Al) 



27r 



47r . 



r 



,,2 / N g din P 
7 ar 

The star radius is the length scale and the dynamical time 
Tdyn = {R^ /GMy^'^ is the time scale. In addition, we ex- 
pand the spheroidal eigenvectors ^ on spherical harmonics 
as (Rieutord 1987) 

1=0 m=-l 



where ¥[^{9, 0) denotes the normalized spherical harmonics. 
We thus obtain the same radial equations coupling uin{r) 
and v1y^(r) as De Boeck et al. (1992), namely 



dr 



> dVm 

dr 



9 



(A2) 



where A = £{£+ 1) and 



7 1 — 



2 



Eliminating leads to the equation for u^^i alone (for clar- 
ity, we drop the subscripts £ and m) 



Au 



, d u 



+ 



du 



dr^ 7 1 — dr 



2 r^(l + r^) 2 A 



7(1 — 7-2)2 



7 Lj2 1 . 



:u = 0. 



After some algebra, this last equation may be written 



A ,/l2 ,d4 ,23 du 
Att -I- u + Br it H — r — 



7 dr 



2 5 du 2'^^'" 2 4^^^^ 
7 dr dr2 dr^ 



dr2 



where 

A = 



- 2A, 



B = 



A. 



"f \ ui / 7 
We now expand u{r) in a power-series of r such as 

oc 

u{r) = "Ycy, (q = £ + l,£ + 3,...), 

9 

where we took into account that u{r) = r^^r oc r^+^ near 
the centre. We thus obtain the following difference equation 
for the Cq coefficients 



B--q-q{q-l) 

1 



+ 



-(g + 2) +2(g + 2)(g + l) 

7 



Cq+2 



-F[A- (g + 4)(g + 3)]c,+4 = 0. 
The convergence of the series implies that 

B--q-q{q-l) = Q, 

7 

which leads to the following expression of the subseismic 
eigenfrequencies 



2 _ 2 



A 



7A-(q-l) (g + f) 



In order to have an easy connection with the eigenfrequen- 
cies of the unapproximated case, we adopt q = £ + 1 + 2k' 
and finally obtain 

2 2 A 

^sub 



7 A - (f + 2k') + 2fc' + 1 + ^) 



with (fc' = 1,2, . 
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A2 The anelastic case 

We start now from Eqs. (|l|) and Q 

^ -> dlnp ^ 



(A3) 



where we used the fact that p is a constant for the homoge- 
neous star. Thus, by applying the same formalism as above, 
we obtain the following anelastic radial equations for u and 

V 



2dv 

u-r — ^ —u, 
ar llI-' 



du 
dr 



(A4) 



and the equation for u alone now reads 

j2.. j2 



Ku — Ar M 1 



'yco^ j dr^ dr^ 
As in the subseismic case, we look for a power-serie so- 

oo 

lution u{r) = ^ ^ Cqr' which leads to the following difference 
equation 

q{q -1)-a(i-—, 
\ "/^ 

+ [A- (q + 2)(Q + l)]c,+2 =0. 

The convergence of the serie requires that 

2 



g(g-l)-A l 



0, 



and we deduce the anelastic eigenfrequencies of the homo- 
geneous star as 

A 



2 _ 2 



7 A-g(g- 1) 
As before, we adopt q — I + 1 + 2k' and obtain 
2 A 



2 



-y A - {£ + 2k'){£ + 2k' + 1)' 
with {k' = 1,2, . . .)• 



where Ma and AAb denote two differential operators. Here 
H is the eigenvector associated with the eigenvalue a^; it 
may read, for instance. 



P 

5P 



These equations are clearly not well adapted for an eigen- 
value problem formulation since Eqs. (|l|),(^) and (^) lead to 
three lines of zeros in the matrix Mb making the system 
not well-conditioned for iterative solvers. 

We therefore prefer to use the oscillation equations of 
Pekeris (1938) who obtained, after judicious substitutions, 
the following reduced system (its equations 12,14 and 15; 
see also Hurley, Roberts & Wright 1966) 



2 dw dg 2dX d* 



^1 



X-^^2^]^^ic^X^gn, + n (B2) 



dr 



d^* 2 d* 



+ - A— -47rG 



(-^+P^)=0, (B3) 



dr^ r dr 



where w = v ■ Cr (should not be confused with the dimen- 
sionless eigenfrequencies uj), X — V • « and $ = ia^p, ijj 
being the gravitational potential perturbation. 

In order to obtain the simplest dimensionless polytropic 
equations, we choose two new length and time scales. As 
length scale, we take R/xi_ whereas To — {A-KGpc)~^^^ is 
our time scale. Here x\ is related to Pc and pc by (see e.g. 
Hansen & Kawaler 1994) 



47rG pi ' 



where Pc and pc respectively denote the central pressure 
and density of the polytrope. With these scales, we have for 
example the following relations 



APPENDIX B: THE EIGENVALUE PROBLEM 
FOR THE POLYTROPE TV = 3 

In this appendix, we formulate the oscillation equations as a 
generalized eigenvalue problem and discuss the 'degeneracy' 
of these equations at the star surface when the subseismic 
or anelastic approximation is applied. 



Bl The complete case 

The governing equations we need to solve are Eqs. (1-4) with 
boundary conditions (6-8). The aim is to derive a generalized 
eigenvalue problem of the form 

XaH = a^MB=., 



2 2 a;iGM 2 R d9 

a = AnGpcA = p;^CJ , g = 47rGpc-^, 

q it^ xi dx 

^ (i^)\.Gpc-^o, q^-m . 

\XlJ 71+1 \dxJ:c = xi 

We note that for the polytrope n — 3, xi ~ 6.89685 and 
q ~ 4.24297 x 10"^ Moreover, we recall that the function 
6{x) is solution of the Lane-Emden equation 

0(0)^1, e{x^) = 0. (B4) 



1 d / ode 

x^ dx 



We thus obtain the dimensionless system (the prime is 
for the derivative d/dx) 



8 B. Dintrans and M. Rieutord 



r + 2-iw 

X I n + 1 



ex' 



+ (1 -7)6»'X - 



■.x' f , 2 

— - \ X — w w 

A \ X 



(B5) 



' w + 



7 



n + 1 

*" + - - ne'e"-^w - e"x = o 



ex + ^, 



A 



with the following set of boundary conditions 
k; = xI/ = at X ^0, 
^ + 1, 



(B6) 



Xl 



-\E' = at X — XI. 



Systems (B5) and (B6) may still be formally written as 



Ma^ = X''MbS., 



(B7) 



where Ma and AAb denote two differential operators and 
A'^ are the real eigenvalues associated with the eigenvectors 
S such as 



w 
X 



B2 The subseismic and anelastic cases 

We have shown in Section 3.1 that the Cowling approxima- 
tion is necessary when the subseismic or anelastic approxi- 
mation are used; thus we now have to deal with 



a^^= N\rer = VX + N'^rCr, (B8) 



p 



c ar 



(B9) 



where we defined x = P' / P- By using the same scales as 
above and eliminating v with Eq. (B9), we obtain the fol- 
lowing dimensionless system for u and p 



A u = X p + \ n 



A \ u + a—u 



Ap. 



(BIO) 



(Bll) 



Here the parameter a differs according to the approximation 
used, namely 

Q — ^ : subseismic, a = — : anelastic, 

7 7 

whereas p = pl^ is related to the projection of x on the 
spherical harmonics as 

At the surface, we have Q(x\) = thus regular singu- 
larities appear in Eqs. ( BIC ) and (Bll). To remove them. 



a solution would consist in multiplying these two equations 
by 6 to obtain the new system 



X^Ou = 6x^p' + [ n 



n + 1 



X^ [Ou + ae'u) = AOp. 



(B12) 



(B13) 



The singularities indeed disappear but are now replaced by 
a surface degeneracy since each equation gives the same so- 
lution u = Q a.t X = x\ (which is however physically correct 
when using the subseismic or anelastic approximation; see 
Section 3.1). This is however making the discretized eigen- 
value problem singular and this should be avoided. 

Following Hurley et al. (1966), we preferred to use the 
Frobenius method. We first develop 9{x) near x = xi &s 

e{x) = e{xi) + {x~ xi)9'{xi) + 0{x - xif, 

and, as 9{xi) — 0, 

e' 1 



for a: — » Xl. 

y X — Xl 

We next expand u{x) and p{x) in series of the form 



(B14) 



u(x) = ^ aq{x - XiY, 

q=0 



p{x) = ^ bq{x - xiy 



q = 



and obtain, using Einstein's notations, 
u = qaq(x — Xl) 



Thus Eq. (Bll) reads now 



A \qaq{x — xi^ + aaq{x ~ xi 



Abqix — XlY 



Equating the lowest power of [x — xi), which is here {x — 
xi)~^ , gives again 

ao = It = 0, 

whereas the next power of {x — xi) gives 

X^ {ai + aai) = Abi A^(l -f- a)u' = Ap. 

We therefore obtain the second boundary condition and 
the system we need to solve is 



2 / , 
X p + 



n+1 

7 



A^ ( It' + a—u I = Ap, 



(B15) 



^ u — Q and A^(l -I- q)ii' = Ap at x = xi. 



which can be rewritten as (B7) with 



